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ABSTRACT 

Polycyclic Aromatic Hydrocarbon (PAH) emission features dominate the mid-infrared spectra of 
star-forming galaxies and can be useful to calibrate star formation rates and diagnose ionized states of 
grains. However, the PAH 3.3 /im feature has not been studied as much as other PAH features since 
it is weaker than others and resides outside of Spitzer capability. In order to detect and calibrate the 
3.3 /im PAH emission and investigate its potential as a star formation rate indicator, we carried out 
an AKARI mission program, AKARI mJy Unbiased Survey of Extragalactic Survey (AMUSES) and 
compare its sample with various literature samples. We obtained 2^5 jam low resolution spectra 
of 20 flux-limited galaxies with mixed SED classes, which yields the detection of the 3.3 fim PAH 
emission from three out of 20 galaxies. For the combined sample of AMUSES and literature samples, 
the 3.3 /im PAH luminosities correlate with the infrared luminosities of star- forming galaxies, albeit 
with a large scatter ( 1.5 dex). The correlation appears to break down at the domain of ultra- luminous 
infrared galaxies (ULIRGs), and the power of the 3.3 /im PAH luminosity as a proxy for the infrared 
luminosity is hampered at log[LpAH3.3/erg/sec~^] > ~42.0. Possible origins for this deviation in the 
correlation are discussed, including contributions from AGN and strongly obscured YSOs, and the 
destruction of PAH molecules in ULIRGs. 
Subject headings: galaxies: star formation 



1. INTRODUCTION 

Over the past decades, infrared (IR) astronomy made 
a tremendous progress largely thanks to various space 
missions, suc h as t he Infra- Red Astronomical Satellite 
(IRAS: Soife r et al.|[l987l). the Infrared Space Observa - 
tory (ISO; K essler et al.lll996l : iGenzel Cesarskyl |2000^ 
and the Spitzer space telescope ([Werner et al.l 12004). 
The most important aspect carried by IR wavelengths is 
probably that IR emission represent s the dust-obscured 
star formation activity of galaxies (jGenzel fc Cesarskvl 
I2OOOI ). The bolometric IR luminosity of galaxies mea- 
sures the dust-obscured star formation within galaxies 
and is less affected by extinction while other shorter 
wavelength star formation rate (SFR) proxies falling on 
optical and ultra-violet (UV) wavelength regimes. How- 
ever using the bolometric IR luminosity as a star forma- 
tion (SF) indicator has two caveats. First, not only newly 
formed massive stars, but also evolved stellar popula- 
tions can heat dust components within galaxies (Calzetti 
[201 li) . Then it is extremely tricky to obtain the whole 
range of IR spectral energy distribution (SED) for high-z 
galaxies to measure bolometric IR luminosities. There- 
fore, many studies have attempted calibrate several IR 
SFR proxy candidates including 8 /im and 24 /im bands 
to bolometric IR luminosities and other SFR indicators 
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in order to establish reliable IR SFR proxies in recent 
years. 

Polycyclic aromatic hydrocarbons (PAHs) have gotten 
enormous attention due to their ubiquity and strong po- 
tential as diagnostics of other properties such as ionized 
states and sizes of grains. PAHs are considered to be 
present in a wide range of objects and environments, 
such as post-AGB stars, planetary nebulae, HII regions, 
reflection nebulae and the diffuse interstel lar medium 
(jPuget et all 119851 : lAllama ndola et all ll989\ The PAH 
features are believed to contribute up t o 10 % o f the^tal 
IR luminosity of star forming; g:alaxies (|Helou et al.ll200ll : 
iPeeters et 311120021 : [Smith et al.ll2007 f). 

Numerous recent studies measure PAH band fluxes and 
equivalent widths (EWs) in order to calibrate these emis- 
sion features as SFR proxies within the Galactic envi- 
ronments and galaxies at higher redshift. These stud- 
ies reveal that there exist differences in PAH EWs and 
L PA H/Lm ratios between local values and hi^h redshift 
ones (|Helou et al.ll200ll : iPeeters et all 120021 : ISmith et all 
|2007[ ). Since PAH band ratios reflect variations in phys- 
ical conditions within environments, such as ionization 
states of dust grains and metallicity (ISniith et al.l 12007 : 
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GalharSHonl), more detailed study on this subject will 
put a better constraint on physical conditions of PAH 
emission sites and calibration of PAH bands as SFR prox- 
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les. 

These previous studies on PAH emission features con- 
centrate on stronger bands, such as 6.2, 7.7, and 11.3 /im. 
On the other hand, studies on 3.3 jam feature have been 
far fewer, due to its relatively weaker strength. There 
still have been efforts to investigate the 3.3 /im PAH 
emission and its property relate d to star formation ac - 
tivity. The first of such studies is lTokunaga et al.l (|l99lf ). 
Compiling inf rared L-band ( 3 -4 jam ) spectra from vari- 
ous sources, Tokunaga et al.l (|l99lf ) categorized the 3.3 
/im PAH emission feature into two types and investi- 
gated if the origin of the emission feature is PAHs. Sev- 
eral studies on energy source of IR emission followed. 
Analyzing 5 7 AG N and one starburst galaxy spectra, 
iClavel et aP (2000) detected the 3.3 /im PAH emission 
from 47 of them an d used its strength for cla ssifying 
the target spectra. iRodriguez^ Ardila ViegasI (L2003) 
detected the 3.3 jam emission in two Seyfert 1 galaxies 
and one quasi- stellar object (QSO) using NASA infrared 
telescope facility (IRTF) SpeX and claimed that these 
active galactic nuclei (AGN) have the 3.3 jam PAH lu- 
minosity l evels simil ar to those of starburst an d LIRGs. 
Ilmanishil (12003') and llmanishi fc Wadal (j2004[ ) observed 
32 Seyfert 2 galaxies and 23 Seyfert 1 galaxies, respec- 
tively, in order to investigate connection between nuclear 
starburst activity and AGN activity using ground-based 
L-band spectroscopy. They detected the 3.3 /im PAH 
emission from 10 out of 23 Seyfert 1 galaxies and 11 
out of 32 Seyfert 2 galaxies respectively, and found that 
starburst activity correlates with nuclear activity regard- 
less of types of Seyfert galaxies. Using ISO/SWS spec- 
tra of a wide variety of sources, Ivan Diedenhoven et aD 
(|2004l ) claimed that the 3.3 /im emission originate mainly 
in neutral and/or negatively charged PAHs in contrast 
to the 6.2 and 7.7 /im emission whi ch are from PAH 
cations. On the other hand, Watabe et al.l p008) and 
lO^i et al.l (|2010f ) confirmed that there is a strong correla- 
tion between nuclear starburst activity and AGN activity 
traced by either X-ray luminosity, or nuclear N-band lu- 
minosity for Seyfert galaxies, while detecting the 3.3 /im 
PAH emission from three out of eight sources and five out 
of 22 sources, respectiv ely. Ex tending these works into 
more luminous regime, llmanishi et al.l (|201ll ) observed 
30 PG QSOs and detected the 3.3 /im PAH emission 
from five QSOs utilizing slit-less spectr oscopic c apability 
of AKARI infrared satellite (Muraka mi et al) 2007) to 
probe global star formation activity by the 3.3 /im PAH 
emission feature.. They confirmed that the correlation 
between nuclear starburst activity and AGN activity are 
intact for PG QSOs. 

There have been other studies utilizing slit-les s spec- 
troscopic capability of AKARI. Imanishi et al.l ()2008D 
obtained AKARI Infrared Camera (IRC; lOnaka et al.l 
l2007l ) spectra of 45 nearby ULIRGs in order to inves- 
tigate the energy source of ULIRGs. They detected 
the 3.3 /im PAH emission from 40 ULIRGs. However, 
they claimed that the obscured starburst activity is not 
the dominant energy source for these ULIRGs, even for 
ULIRGs which are classified as non- Seyfert optically. 
llmanishi et al.l (|2010l ) also obtained IRC spectroscopy of 
64 LIRGs and 54 ULIRGs and detected the 3.3 /im PAH 
emission from the majority of the observed targets. They 
found that the 3.3 /im PAH gives a good estimate for 
SFR as Bra {Xrest = 4.05 /im) does. For more of higher 



redshift objects. fSajina et al.l (|2009f ) detected the 3.3 /im 
PAH emission for four out of 11 z ~ 2 ULIRGs. 

Understanding how the 3.3 jam PAH emission is re- 
lated to star formation activities has a great importance, 
since the 3.3 /im PAH feature is likely to be the only 
dust emission feature at high redshift {z > 4.5) easily 
accessible with future space IR telescope missions, such 
as James Webb Space Telescope (JWST) and SPICA. 
Given the importance of understanding 2.5-5 fim mid- 
infrared (MIR) emissions from extragalactic sources, we 
carried out a study of low redshift objects in the 2.5- 
5 fim window, AKARI mJy Unbiased Survey of Extra- 
galactic Sources (AMUSES) as one of AKARI mission 
projects (MPs). The main scientific goal of AMUSES is 
to construct a continuous spectral library over the wave- 
length window between 2.5 and 40 fim for a subsample of 
pmJy Unbiased Spitzer Extragalactic Survey (5MUSES, 
IWu et aLll2010f ) by combining spectra from AKARI and 
Spitzer. 

In ^ we present sample selection, data acquisition and 
reduction. Then we present and discuss fitting method- 
ology as well as stacking analysis in ^ In 21 we present 
both reduced individual spectra and the stacked spectra 
of the sample. We also present and discuss fitting results 
in 21 Implication of our results are presented in 21 We 
summarize our study in 21 

Throughout this paper, we assume that the universe 
is flat with flM= 0.3 and Q\=0.7, and i7o=70km s~^ 
Mpc-^ fe.g.. llm et allll997[ ). 

2. DATA 

2.1. Sample Selection 

The AMUSES sample is drawn from 5MUSES. We 
choose 5MUSES as our parent sample, because it con- 
tains a statistically unbiased library of MIR spectra of IR 
sources that covers the gap between the bright, nearby IR 
galaxies and the much fainter, more distant IR sources 
being studied by the Spizter/ AKARI . 5MUSES, one 
of the Spitzer Legacy surveys, performed a MIR spec- 
troscopic observation of extragalactic sources brighter 
than b mJy at 2Aii m in the Spi tzer Fi rst Look Survey 
(FLS ; fCondon et a l. 2003; Fad da et al.i r2004: Y an et al 



2004; 



2006 



Appleton et al 2004: Choi et al.l [2006: Shim et al 



2007f ) field and four subfields of Spi tzer Wide-area 



Infrared Extragalactic survey (SWIRE; "Lonsdale e t al.l 
2004) with the Infrared Spectrograph (IRS) onboard the 
Spitzer space telescope. The main scientific goal of 
5MUSES is to provide an unbiased library of infrared 
spectra from 5 to 40 /im of sources which have not been 
sought after in previous studies. 

Since the main objective of AMUSES is to detect the 
3.3 pm PAH feature, we select 5MUSES galaxies that are 
brighter than 1 mJy at 3.6 jam. Based on the Spitzer 
IRAC 3.6 /im data for the 5MUSES sample, we find that 
60 of the 330 5MUSES galaxies satisfy the flux cut of 
ImJy at 3.6 jam. We limit the redshift range to z < 0.5 
in order to achieve S/N > 5 for detection of the 3.3 /im 
emission feature. With this additional cut, the sample 
size is reduced to 50 galaxies. In addition to this base 
sample, we add 10 5MUSES galaxies with their 3.6 /im 
flux brighter than 0.7 mJy whose redshifts could not be 
determined by optical spectroscopy. These are added as 
secondary targets in case that the scheduling constraint 
makes it impossible to observe the main targets. With 
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these additional galaxies, the original sample consists of 
60 galaxies. 

Ultimately, a total of 44 galaxies among these 60 galax- 
ies are approved for the final program, after resolving 
visibility conflicts with targets from other approved pro- 
grams. However, only 20 target galaxies were observed 
until the mechanical failure of the AKARI cooling sys- 
tem. The termination of the AKARI scientific mission 
is announced on 2011 June 20, thus, these 20 galaxies 
are the final sample for AMUSES. F or our fina l sam- 
ple, we adopted their SED classes from lWu et al.l (pOlO). 
They classified objects using the 6.2 /im PAH equivalent 
widths (EWs): sources with EW > 0.5 /im are starburst 
(SB)-dominated; sources with 0.2 < EW < 0.5 jam are 
AGN-SB composite and sources with EWs < 0.2 fim are 
AGN-dominated. There are 12 AGNs and seven star- 
burst galaxies within the final sample. One galaxy has a 
composite SED. The final sample also includes two AGNs 
at z > 0.7 whos e redshifts were d etermined after our sam- 
ple selection by'Wu et al.' ("2010). The basic properties of 
the sample are listed in Table [TJ while the redshift versus 
the 3.6 jiim flux is presented in Figure [H 

2.2. Observation 

The observation was carried out with the spectroscopy 
mode of IRC on AKARI. The NIR grism (NG) mode was 
adopted for most of the ob servations providin g the spec- 
tral resolution of R - 120 (jOnaka et al.l l200f). For eight 
of the targets, NIR prism (NP) mode was used which 
provide the spectral resolution of R ~ 20. The observa- 
tions were carried out using 1' x 1' slit aperture to avoid 
confusion with surrounding sources. The pixel scale of 
NIR camera is 1.45'^ and the full width half maxima of 
the point spread function is 3.2 pixels which corresponds 
to 4.64'^ All target galaxies have sizes small enough to 
fit in the slit aperture, although size variation is not in- 
significant. We present the R-band images of the target 
galaxies along with the ID reduced spectra and the 2D 
spectrum images in Figure [2j For several galaxies with- 
out ancillary R-band images, Spitzer IRAC 3.6 /im band 
images are presented. We discuss the extraction width 
and the target galaxy sizes in section [231 The actual ex- 
posure times run between six to seven minutes for each 
pointing observation. Originally, three pointing observa- 
tions were planned for all the targets, but a half of the 
targets presented in Table 1 had less than three pointing 
observations due to the termination of the AKARI sci- 
ence mission. In total, observation of 51 pointings was 
carried out for 20 target galaxies. Table [1] also lists how 
many pointings are acquired for each target. 

2.3. Data Reduction 

Data reduction was performed with the IRC spec- 
troscopy pipeline.^ The IRC spectroscopy pipeline sub- 
tracts scaled super-dark frames, applies linearity cor- 
rection while masking saturation and monochromatic 
flat-fielding, subtracts background from individual sub- 
frames, registers and stacks images before subtracting 
background again from stacked images and extracting 
2D spectra. Then wavelength calibration, flat color-term 
correction, spectral tilt correction, and spectral response 
calibration were applied on these extracted 2D spectra. 
Additional cosmic ray removal, stacking multi-pointing 
exposures, and one-sigma clipping during stacking were 



executed individually after processing the data through 
the IRC spectroscopy pipeline. 

Final one-dimensional (ID) spectra were extracted 
from two-dimensional (2D) spectral images which are 
binned by three pixels along the wavelength direction. 
The ID spectra were extracted from the 2D spectrum 
images with various widths in spatial direction. First, we 
extracted the ID spectra with the width of 5 pixels which 
corresponds to 7.3" . In physical scales, this ranges from 
1.8 kpc for the nearest target galaxy to 56.3 kpc for the 
farthest target galaxy. We present the extraction boxes 
with 5 pixel width over the R-, or IRAC 3.6 jam band 
images in Figure [2j As indicated in these figures, these 
extraction regions cover entire galaxies in most cases. We 
tried to change the extraction width in order to check if a 
larger extraction width unveils any hidden PAH 3.3 /im 
emission feature from them, or it changes the extracted 
PHA flux more than the measurement errors. Specif- 
ically, we also applied extraction boxes with 10 pixels 
and 15 pixels. However, we did not notice any signif- 
icant change of the extracted spectra over the range of 
extraction width for our entire sample. Therefore, we use 
only the spectra extracted with 5 pixel extraction width 
for our analysis. 

For more details regarding data reduction process, 
please refer .Ohv ama et al. (2007). 

3. ANALYSIS 
3.1. Spectra of the sample and their analysis 

We present the reduced spectra and fitting of emis- 
sion features of 20 AMUSES galaxies in Figure O For 
FBQS J0216-0444 and SDSS J160128.54+544521.3, the 
presented spectra are based on Spitzer IRS spectra of 
5MUSES due to high redshifts instead of AKARI IRC 
spectra. We also include insets showing fitting results of 
the 3.3 /im PAH features for individual galaxies detected 
with the 3.3 jam PAH emission. 

In order to detect the 3.3 jam PAH emission feature 
and measure its strength, we fit the spec tra a round 3.3 
/im usmg the Drude profile (|Li fe Drainell200ll) : 



brlr 



{\/K-\rl\Y+l^ 



(1) 



where is the central wavelength of the feature, 7^ is 
the fractional FWHM, and hr is the central intensity. 

Fitting and subtracting continuum from spectra also 
can be a subject of debate. The most difficult issue is 
how to avoid absorption features which are abundant 
throughout MIR wavelength regimes. These absorption 
features include CO2 ice absorption at 4.27 /im, ^^C02 
ice absorption at 4.38 /im, and CO ice absorption fea- 
ture at 4.67 /im. The vicinity of the 3.3 jim PAH feature 
is not an exception with water ice absorption feature at 
3.05 /im and carbonaceous dust absorption feature at 3.4 
/im. Therefore, we decide to apply a spline fit for con- 
tinuum to the wavelength range between 3.0 jim and 3.6 
pm while masking out the 3.3 /im emission feature and 
avoiding any obvious absorption feature by visual inspec- 
tion. Generally, we used the wavelength ranges between 
3.0 /im and 3.2 /im, and between 3.4 /im and 3.6 /im 

^ http:/ /www. ir.isas.jaxa.jp/ASTRO- 

F/Observation /DataReduction /IRC / 
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for the continuurQ fit. If the 3.3 /im feature looks to be 
shifted upon visual inspection (perhaps due to an un- 
known systematic error in wavelength calibration), we 
adjust the fitting ranges. However, there is no 3.3 /im 
feature detected with a shifted peak. Also, if there is a 
clear absorption feature residing within the ranges upon 
visual inspection, we adjust the fitting range further out 
from 3.3 /im. 

We use an IDL procedure, MPFIT, to apply the Drude 
fit onto the 3.3 jam PAH feature after subtracting contin- 
uum from each individual spectrum. We do not constrain 
peaks and widths of Drude fit unless the IDL procedure 
cannot fit any obvious emission feature. We confirm the 
detection of the 3.3 jam PAH emission feature by visual 
inspection of the fitting, while considering S/N over the 
continuum. We measure line fluxes and EWs of the 3.3 
jiim PAH emission based on the outputs of the fitting 
results of Drude fits and continuum fits for the sample 
galaxies. The flux errors are determined based on the 
Drude fit and its errors measured by MPFIT. These val- 
ues are given in Table [2l 

When no 3.3 jam PAH emission is detected, we mea- 
sure upper limits. In order to calculate the upper limits, 
first, we measure the standard deviation of data points 
within a wavelength range at which the 3.3 /im should 
be located. This range is decided by the average width 
obtained from the other galaxies detected with the 3.3 
PAH jam emission. The average width is 17 nm. Then 
we assume a Drude profile which has a peak three times 
bigger than the standard deviation, and has a width of 17 
nm which is the average width of the 3.3 jam PAH emis- 
sion detected within our sample. We take the sums un- 
der these Drude profiles as upper limits for non-detection 
targets. 

3.2. Stacking analysis 

In order to construct representative spectra of each 
SED class while recovering missing 3.3 jim features due 
to low S/N, we apply stacking analysis. We constructed 
blind-stacked spectra of the entire sample, the galax- 
ies with the 3.3 jim PAH feature detection, the galaxies 
without the 3.3 jim PAH feature, the galaxies with AGN 
SEDs, and the galaxies with starburst SEDs. 

To construct the stacked spectra, we simply repro- 
duced spectra of individual sample galaxies onto a wave- 
length grid set to have even step through the wavelength 
range between 2.5 jam and 8 /im, then constructed the 
spectra for each group by simply averaging spectra of 
group galaxies without applying any clipping or normal- 
ization. Spitzer IRS spectra provide data for any wave- 
length range required to construct the new spectra. We 
present the stacked spectra of each galaxy group in Fig- 
ure [3l 

We concern that the stacked spectra are dominated by 
a few galaxies with high flux. Therefore we also stacked 
normalized individual spectra for each group and com- 
pared the spectra to the stacked spectra produced with- 
out normalization. To do this, we normalized individual 
spectra by fluxes between 3.5 jim and 4.0 /im, then ob- 
tained averages of them. However, overall, these stacked 
spectra do not show any significant departures from the 
ones in Figure [3l 

4. RESULTS 



4.1. Detection and strength of the 3.3 jim PAH 

emission 

We detect the 3.3 jim PAH feature from three galax- 
ies out of 20 target galaxies. Although the 3.3 jim PAH 
features are generally pretty weak across the sample, a 
couple galaxies of starburst SED class, namely 2MASX 
J16182316+5527217 and 2MASX J16205879+5425127 
show very clear and strong 3.3 jim PAH features. While 
these two galaxies have starburst SEDs, the remaining 
one galaxy, 2MASX J10542172+5823445 has an AGN 
SED. Therefore, the detection rate for the starburst SED 
class sample is 29% (2/7), while the detection rate for the 
AGN SED class sample is 8% (1/12). However, since the 
desired S/N is not achieved for a large portion of the sam- 
ple, these detection rates should be considered as lower 
limits. 

Upper limits varies significantly. However, generally it 
is smaller for the starburst SED galaxies than the AGN 
SED galaxies. While the five SB galaxies without the 
3.3 jim PAH emission detection have their average upper 
limits of 5.97 X 10^^ erg s~^, their AGN counterparts have 
detection limits of 1.35x10^^ erg s~^. 

The 3.3 jim PAH luminosities range from 7.11 
xlO'^^ erg s-^ of an AGN at z = 0.205 (2MASX 
J10542172+5823445) to 7.14 xlO^^^ erg s"^ of a star- 
burst galaxy at z =0.063 (2MASX J16205879+5425127). 
The two starburst galaxies with 3.3 /im PAH detection 
are at z ~ 0.1 and their average luminosity is 9.46 xlO^^ 
erg s~-^. 

4.2. Stacked spectra 

The stacked spectrum of the entire sample shows a 
distinctive emission feature of the 3.3 jim PAH emission. 
The EW of the stacked spectrum is 11 nm, while the 
spectrum stacked with the galaxies with starburst SEDs 
has the EW of 32 nm. On the other hand, the spectrum 
stacked with the non-detection galaxies does not reveal 
the 3.3 jim PAH emission feature. The 3.3 jim PAH 
emission does not show up for the stacked AGN spec- 
trum, either. It is quite obvious that AGN host galaxies 
do not have strong PAH emission overall. 

5. DISCUSSION 

5.1. The Correlation between Lpah3.3 ci'rid Ljr 

We present a plot comparing luminosity of the 3.3 jim 
PAH emission (Lpahs.s) with infrared luminosity, Lxr 
(Figure Hj). Due to lack of the detected sources within the 
AMUSES sample, it is not meaningful to derive the cor- 
relation between Lpah3.3 and Ljr only for the AMUSES 
sample. Therefore we look for the correlation for the 
combined sample of AMUSES and literature samples. 
Big filled circles represent samples of the AMUSES ob- 
jects, while small symbols represent san i ples from liter- 
atures ( Rodriguez- Ardila fc Viegasl'2003': Imanishi e t al. l 
2008; Saiina et al. 2009; Imanishi et al. 2010; Lee e tHI 
2012). Note that Ljr come from Wu et al. (2010). In or- 
der to estimate Lir of the sample of 5MUSES, Wu e t al.l 
(|2010l ) utilize mid-IR spectra from IRS along with the 
IRAC and MIPS photometry. First, they construct an IR 
template library based on these data. This library covers 
a wide range of galaxies from normal star-forming galax- 
ies to ULIRGs to quasars. Then using synthetic IRAC 
photometry drawn from the IRS spectra and MIPS pho- 
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tometry, they compare them with the corresponding syn- 
thetic photometry from the SED templates and estimates 
LiR. On the other hand, the other samples utilize IRAS 
photometry to calculate Lir. Mor e specifically, they 
used t he formula given in Table 1 of ISanders fc Mirabell 
([19961 ). 

The figure shows a broad correlation between Lpah3.3 
and I/iR for the combined sample of our data points and 
the data from the literature, though the scatter in the 
relation is large. In order to obtain a linear correlation 
between Lpah3.3 and Lir, we carry out fitting to the 
combined sample while switching abscissa and ordinate 
between Lpah3.3 and Ljr, since the slope of the correla- 
tion changes significantly (by 15 %) depending on which 
quantity is chose as a dependent. Then we take the aver- 
age fit as our fitting result. The linear fit to the combined 
sample is as in Eq. (j2j); 

log(LiR) = (1.16 ± 0.30) X log(LpAH3.3) 

-(3.11 ±0.34), (2) 

where Ljr and i^pAH3.3 are in the unit of erg sec~^ with 
the correlation coefficient of 0.70. This fit is shown in 
Figure H] with a solid line. The undetected sources are 
not included for this fitting. In general, the data points 
of AMUSES sample show higher Lpah3.3 than the liter- 
ature sample. Considering the fact that we may have 
a rather shallow detection limit and an unbiased dis- 
tribution for AMUSES sample, the detected sources of 
AMUSES sample are more luminous objects in the 3.3. 
jam PAH emission feature and occupy the upper envelope 
in terms of I/pah3.3 for the correlation between Lpah3.3 
and I/IR of the literature samples. 

On the other hand, Ljr can contain non-star-forming 
contribution. Therefore, we consider the correlation be- 
tween LpAH3.3 and Ljr only for sources with starburst 
SEDs, or Hll-like SEDs. The fit is given as in Eq. (|3]); 

log(LiR) = (1.37 ± 0.17) X log(LpAH3.3) 

-(12.18 ±0.75), (3) 

where the units for Lpah3.3 and Ljr are same with Eq. 
(|2]). The correlation coefficient for the fit is 0.69. We use 
the same method which we use for Eq. (|2]) for this fit. 
This fit is shown within Figure H] by a dotted line. 

5.2. The PAH 3.3 iim emission as a star formation 
indicator 

In order to to be an effective star formation indicator, 
^PAH3.3 must have a good correlation with Ljr. Within 
the detected sources from AMUSES sample and the liter- 
ature samples, especially with the sources with SB SEDs, 
^PAH3.3 and LiR have a correlation between them, but 
with a large scatter (Figure H]). 

This can be also confirmed by previous studies on other 
PAH emission features. Although systematic calibration 
of any particular PAH emission feature to SFR has been 
carried out less extensively than other SFR indicators, 
there have been a few studies to show strong correla- 
tions between the streng:ths of PAH e mission features 
and IR luminosity. iGenzel et al.l (|l998 f) utilizes the the 
strength of the 7.7 /im PAH emission for diagnosing en- 
ergy sources for ULIRGs based on 15 spectra from ISO. 
Also usi ng ISO spectra of ga lactic and extragalactic star 
sources, iPeeters et al.l (|2004l ) shows that the luminosity 



of the 6.2 jam P AH emission (I/pAHfi .2) correlates very 
well with Lfir- iBrandl et al.l (|2006f ) calibrate the 6.2 
/im PAH emission to Lir using IRS spectra of 22 star- 
burst galaxies. Regar ding more of higher redshift ob- 
jects. Pope et al. (2008) found that all the PAH bands 
which they detected from 13 sub-mm galaxies (SMGs) at 
various redshifts between 1.0 and 2.5, namely the 6.2, 7.7, 
and 1 1.3 /im PAH eni ission features, correlate well with 
LiR. IShi et al] (|2009l ) also found that the 6.2 /im and 
7.7 /im PAH emission features correlate well with SFR 
based on Ljr and use them t o measure the cosm ic star 
formation history of QSOs. iFiolet et al.l (|2010f ) shows 
that the 6.2 /im and 7.7 /im PAH luminosities of their 16 
ULIRGs at z ^ 2 follow the correlation of the samples 
of PoDe et al. (2008), Menendez-Delmestre et al. (20Q9|i, 
and Sh i et a l. (2009). On top of these studies, Wu et al.l 
(j2010f ) found strong correlations between I/pah6.2 and 
LiR for SB, composite, and AGN sources based on 
5MUSES, the sourc e of the parent s ample for this study 
(Al, A2, and A3 of IWu et al.l (j20Toh V 

Whether our I/pah3.3-^ir correlation persists, espe- 
cially at high LiR remains as a subject of further in- 
vestigation. In Figure [5l we plot L t r aga ins t I/pah6.2- 
We a dopt Lpah6.2 from IWu et al.l (l2010h. IWu et al.l 
(j2010f ) utilize the PAHFIT software (Smith et alJ|2007D 
and measured Lpah6.2 from the fitting result based on it. 
A linear fit to the detections of AMUSES sample gives 
the correlation between Lpah6.2 and Ljr as in Eq. (|4]); 

log(LiR) = (0.95 ± 0.03) X log(LpAH6.2) 

+(4.21 ±0.19), (4) 

where the units for Lpah6.2 and Ljr are same with Eq. 
(j2j) with the correlation coefficient of 0.75. However, it is 
noticeable that the ULIRGs deviate from the correlation. 
They have either higher Ljr, or lower Lpah6.2 than lower 
^IR objects suggest. It looks plausible that ULIRGs at 
LiR > 10^^ M0 may deviate from correlations between 
PAH luminosity and Lir, since these ULIRGs also may 
deviate from the correlation between I/p ah3.3 and Ltr 
(Figure |4|) . Similar trends a re rep orted by 'Peeters et al.l 
^004) and Imanis hi et al.l (120071 ). In order to check if 
this is true, we obtain a linear fit to the detections with 
SB SEDs which are not ULIRGs. This fit gives the cor- 
relation between Lpah3.3 and Ljr as in Eq. (|5j); 

log(LiR) = (0.82 ± 0.13) X log(LpAH3.3) 

+(10.58 ±0.86), (5) 

where the units for Lpah3.3 and Ljr are same with Eq. 
(j2j). The correlation coefficient for the fit is 0.71 and the 
fit is shown in Figure |4] as a dashed line. 

It is noticeable that ULIRGs generally have smaller 
^PAH3.3 for a given Ljr regardless of SED class es based 
on th is fit. Similar trends are found by Cla vel et al.l 
(120001) for the 7.7 /im PAH emission feature and 
IPeeters et aP (j2004f ) for the 6.2 /im PAH emission fea- 
ture. It is not plausible that the offset for ULIRGs are 
due to any systematic difference in measurement meth- 
ods for LpAH3.3, or Ljr. For example, it is unlikely if 
LiR based on IRAS four band photome try overest i mates 
true LiR, or Ljr which we adopted from lWu et al.l (|2010l ) 
underestimates. 

There can be two possible explanations why ULIRGs 
show the trend. First, ULIRGs may have larger non- 
star-forming contribution to their Ljr regardless of their 
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SED classes based on their optical spectra. There have 
been studies claiming tha t AGN contribution to Ltr is 
not negligible for ULIRGs (jNetzerll2QQ9l: IShao et al.ljgQaQl . 
and references therein). For instance, iLutz et aTl (| 19981 ) 
show that, while only 15% of ULIRGs at luminosities 
below 2 X 10^^ L0 are attributed to AGN, this frac- 
tion increase s to 5 0% at higher luminosity. Figure 5 of 
iNetzer et al.l (|2QQ7 ) show the correlations between nu- 
clear activity and two IR luminosities; FIR(60/im) lumi- 
nosity (Leo) and the 7.7 jam PAH luminosity. Interest- 
ingly, the correlation between the nuclear activity probed 
by the A5100A monochromatic luminosity and the 7.7 
/im PAH luminosity breaks at about Leo ^ 10^^ L©, 
while Leo correlates well with the entire range of A5100A 
monochromatic luminosity. The difference between these 
two correlations support that AGN activity may con- 
tribute to LiR more than it has been considered to be. 
Also, considering the fact that there is n o sign of depar- 
ture for ULIRG sample of llmanishi et al.l (jJOOS) from the 
correlation between Lpah3.3 and Lpah6.2 (Figure [6j), it is 
more likely to be Lm which causes the break. Recently, 
lYamada et al.l (|2Q12l ) found that the correlation between 
^PAHS.s/^iR and Ljr has a break around Ljr ~ 10^^ Lq 
from another MP of AKARI, MSAGN. There is another 
candidate to contribute to Ljr for ULIRGs other than 
AGN activity: embedded young stellar obje cts (YSOs) 
(jPeeters et al.ll2QQ4l : ISpoon et al.ll2QQ2l I2QQ7D . Although 
embedded YSOs are from star- forming regions, they 
should be considered as a non-typical star- forming; con- 
tribution to Ljr. Based on ISO observation. ISpoon et al.l 
(2002) found that most ULIRGs have the 6 - 8 //m ice 
absorption feature and linked it to strongly obscured be- 
ginning of star formation. 

Second, ULIRGs suppress PAH emissions, not only the 
3.3 /im PAH emission, but the PAH emission features in 
general. As shown in Figure|4]and FigurejSl both Lpah3.3 
and LpAH6.2 are lower for ULIRGs for a given Ljr. It 
is believed that strong UV radiation, such as one from 
AGN, can destruct PAH molecules and reduce the pop- 
ulation of PAHs. In fact, Imanishi et al.l (|2008') claimed 
that the ratio between the 3.3 fim PAH luminosity and 
infrared luminosity for ULIRGs are less than what is ex- 
pected from obscured starbursts and are attributed to 
the deficiency of PAH within ULIRGs. However, it is 
unclear why AGN at Ljji < 10^^ L© does not reduce the 
PAH luminosity as much. Therefore, it is plausible that 
most ULIRGs have a large contribution from AGN to 
their Ljr. A mult i- wavelength analysis of the ULIRGs, 
including the mi d-IR component from hot dust, should 
settle this issue (|Nardini et al.ll2009L 120101 : iSaiina et al] 
[2012). 

However, it is still possible that there is no such depar- 
ture from the correlations for ULIRGs. Since there are 
few sources detected with the PAH emission features at 
the lower end, the correlations between the PAH lumi- 
nosities and LiR can be skewed by the few objects with 
large non-star- forming contribution to their Ljr. In or- 
der to see if it is the case, it is required to detect more 
lower luminous objects in terms of PAH luminosity. 

On the other hand, as shown in Figure H] and Figure |5l 
the 3.3 /im and 6.2 /im emission features behave similarly 
against Lir. In order to see if there is any difference be- 



tween them, we plot Lpahs.s against Lpah6.2 and present 
it in Figure |6l Although it looks to have large scatters, 
the 3.3 /im PAH emission has a strong correlation with 
the 6.2 /im PAH emission. A linear fit to the data points 
of all samples without the upper limits gives the correla- 
tion between Lpah3.3 and Lpah6.2 as in Eq. (|6]); 

log(LpAH6.2) = (0.83 ± 0.06) X log(LpAH3.3) 

+(7.88 ±0.41), (6) 

where the units for Lpahs.s and i^pAH6.2 are erg sec~^ 
with the correlation coefficient of 0.83. This fit is shown 
within Figure [6] by a solid line. Interestingly, ULIRGs do 
not deviate from the overall trend. These ULIRGs devi- 
ate rather significantly from the correlation between Ljr 
and PAH luminosities (Figure |4] and Figure |5j). There- 
fore, it assures that the 3.3 /im PAH can represent what 
the 6.2 /im represents, if needed. However, it is notice- 
able that the correlation between Lpahs.s and Lpah6.2 
is far from unity. This may be attributed to the ori- 
gin of PAH emission f eature s. As mentioned before, 
Ivan Diedenhoven et al.l (|2004 f) compare variations of the 
profiles of various PAH emission features and show that 
while the 6.2 and 7.7 /im emission originate from PAH 
cations, the 3.3 /im emission feature is attributed to neu- 
tral and/or negatively charged PAHs, which is suggested 
by several labo ratory studies jS zczepanski fc Valal 119931 : 
Langhoff 1996; i Hudgins Allam andola 199^^ 

Therefore, although the 3.3 /im PAH emission feature 
shows that it has the potential to be a SFR proxy like 
other PAH emission features, it heavily relies on under- 
standing the physical conditions of PAH molecules and 
their emission features. It also requires to understand 
the deviation of ULIRGs and its origin lest the power of 
^PAH3.3 as a proxy is hampered by this deviation. 

6. SUMMARY 

We describe the details of a mission project of AKARI, 
AMUSES and its results. With AMUSES, we have in- 
vestigated the 3.3 /im PAH feature and its correlation 
with other PAH features as well as other properties of 
the sample galaxies. The summary of our results are: 

1. We detected the 3.3 /im PAH emission from three 
out of 20 flux limited sample galaxies and the de- 
tection rate for our sample is 15%. 

2. For the combined sample of AMUSES and the lit- 
erature samples, Lpah3.3 correlates with Ljr, while 
^PAH6.2 also has a similar correlation with Ljr. 

3. ULIRGs at Lm > IO^^Lq may deviate from the 
correlation between Lpah3.3 and Lir of the lower 
luminous objects due to non- star- forming contri- 
bution on LiR, such as AGN activity and heavily 
obscured YSOs, or destruction of PAH molecules 
by AGN activity. 
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TABLE 1 

Basic Properties of Sample Galaxies 



Galaxy 


RA 


Dec 


redshift 


(mjy) 


SED class 


number of pointings 
NG NP 


FBQS J0216-0444 


2:16:40.72 


-4:44:05.1 


0.87 


2.11 


AGN 


3 


1 


2MASX J02165778-0324592 


2:16:57.77 


-3:24:59.7 


0.137 


2.53 


AGN 


1 





2MASX J02191605-0557269 


2:19:16.05 


-5:57:26.9 


0.103 


1.04 


AGN 


Q 
O 


i 


2MASX J02193906-0511336 


2:19:30.08 


-5:11:33.8 


0.151 


1.93 


AGN 


1 
i 


n 
U 


2MASX J02195305-0518236 


2:19:53.04 


-5:18:24.1 


0.072 


1.40 


SB 


1 
i 


U 


2MASX J10520659+5809476 


10:52:06.56 


+58:09:47.1 


0.1172 


1.34 


SB 


1 
1 


n 
u 


2MASX J10542172+5823445 


10:54:21.65 


+58:23:44.6 


0.2045 


2.58 


AGN 


2 





VII ZW353 


10:57:05.43 


+58:04:37.4 


0.1404 


2.35 


AGN 


3 


1 


SDSS J105959. 93+574848.1 


10:59:59.95 


+57:48:48.1 


0.4530 


1.20 


AGN 


3 





SDSS J155936.13+544203.8 


15:59:36.12 


+54:42:03.7 


0.3077 


1.09 


AGN 


3 


1 


SDSS J160128.54+544521.3 


16:01:28.52 


+54:45:21.3 


0.7278 


1.13 


AGN 


3 


1 


2MASX J16044063+5534089 


16:04:40.64 


+55:34:09.2 


0.078 


1.23 


SB 


2 





2MASX J16130186+5521231 


16:13:01.82 


+55:21:23.0 


0.012 


2.07 


SB 


1 





2MASX J16144902+5545120 


16:14:45.92 


+55:45:12.9 


0.064 


1.55 


AGN 


2 





SBS1614+546 


16:15:21.78 


+54:31:48.3 


0.474 


1.18 


AGN 


1 





SDSS J161445.94+542554.4 


16:16:45.92 


+54:25:54.4 


0.223 


1.06 


AGN 


1 





2MASX J16165997+5600276 


16:16:59.95 


+56:00:27.2 


0.063 


2.35 


SB 


3 


1 


2MASX J16181934+5418587 


16:18:19.31 


+54:18:59.0 


0.083 


1.50 


Composite 


2 





2MASX J16182316+5527217 


16:18:23.11 


+55:27:21.4 


0.084 


1.71 


SB 


3 


1 


2MASX J16205879+5425127 


16:20:58.82 


+54:25:13.1 


0.082 


1.71 


SB 


3 


1 



Note. f Spizter IRAC 1 band fluxes from 5MUSES LWu et al.ll2QrQV 



TABLE 2 

PAH MEASUREMENT OF THE SAMPLE GALAXIES 



Galaxy 


/3.3 


EW3.3(nm) 


J6.2 


J7.7 


J 11.3 


log LIR(M/M0) 


FBQS J0216-0444 


<589.9 


<10.4 


2.407 




15.23 


12.70 ± 0.01 


2MASX J02165778-0324592 


<4320.3 


<38.2 


1.108 


6.167 


10.04 


10.90 ±0.03 


2MASX J02191605-0557269 


<1853.6 


<9.3 


14.57 


59.89 


33.77 


10.71 ± 0.05 


2MASX J02193906-0511336 


<3268.3 


<49.8 


44.66 


230.2 


72.83 


11.38 ±0.06 


2MASX J02195305-0518236 


<3504.7 


<17.1 


292.7 


931.5 


253.5 


10.93 ±0.03 


2MASX J10520659+5809476 


<5290.1 


<3.2 


231.3 


821.1 


228.6 


11.34 ± 0.03 


2MASX J10542172+5823445 


5882.3+2083.9 


115.2+40.8 


50.21 


222.8 


65.59 


11.43 ± 0.03 


VII ZW353 


<1703.0 


<4.9 


40.10 


220.6 


92.47 


11.18 ± 0.03 


SDSS J105959.93+574848.1 


<2873.6 


<32.3 


5.993 


20.28 


11.79 


11.83 ± 0.02 


SDSS J155936.13+544203.8 


<2795.4 


<56.7 


6.254 


13.40 


22.42 


11.32 ± 0.06 


SDSS J160128.54+544521.3 


<794.4 


<18.5 


7.008 


17.24 


41.84 


12.47 ± 0.01 


2MASX J16044063+5534089 


< 2608.4 


< 4.5 


152.5 


559.9 


173.7 


11.10 ± 0.04 


2MASX J16130186+5521231 


<6279.6 


<10.5 


291.1 


1045 


338.4 


9.47 ± 0.05 


2MASX J16144902+5545120 


<3712.7 


<7.2 


46.75 


103.7 


35.10 


10.26 ± 0.03 


SBS1614+546 


<7860.7 


<125.3 


5.825 


16.20 


12.97 


11.47 ± 0.08 


SDSS J161445.94+542554.4 


<3484.4 


<43.9 


33.79 


88.38 


45.77 


11.26 ± 0.02 


2MASX J16165997+5600276 


< 2574.6 


< 2.2 


165.0 


597.8 


185.5 


10.66 ± 0.02 


2MASX J16181934+5418587 


<3975.1 


<4.8 


173.4 


685.4 


193.6 


11.14 ± 0.04 


2MASX J16182316+5527217 


6734.2+2053.4 


170.1+51.9 


435.4 


1456 


398.5 


11.13 ± 0.03 


2MASX J16205879+5425127 


4290.7+1807.3 


67.3+28.3 


297.6 


954.9 


259.1 


11.11 ± 0.03 


Note. — f Line fluxes of the 3.3 fim PAH feature 


with units of 10 ""^^erg 

- 1 _ 1 _o 


s~"'^cm~ 


ij: Line fluxes of the other 



PAH feature from lWu et al.l (l2Qia ) with units of 10~^^ erg s~^cm~^. 
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Fig. 1. — Redshift distribution of the sample of galaxies of AMUSES plotted against 3.6 /im flux. Symbols represent the SED classiflcation 
of target galaxies based on their MIR colors by 5MUSES collaboration. Filled circles represent the AGN-type SED, while asterisks represent 
the starburst-type SED. An open diamond represents the composite SED target. AGN-type targets are more evenly distributed across the 
redshift range, while starburst-type targets are more clustered at lower redshifts. 
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Fig. 2. — Continued. 
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Fig. 2. — The spectra of AMUSES based on IRC and IRS spectra along with ancillary images. All spectra are de-redshifted based 
on their redshifts, therefore x-axis is rest-frame wavelength. For two galaxies at high redshifts, FBQS J02191605-0557269 and SDSS 
J160128. 54+544521. 3, the spectra presented are based on IRS spectra. For the rest of the sample, the spectra are based on IRC spectra. 
Each panel includes the name of a target galaxy and its SED class as well. For galaxies detected with the 3.3 fim PAH emission feature, 
we mark the locations of the feature by upward arrows and also show fitting results. Fitting results are shown with blue lines and insets. 
Overplotted blue lines are sums of continuum fits and fitted Drude profiles. Insets show the 3.3 fim PAH emission feature after the fitted 
continuum subtracted from the spectrum. Then overplotted red lines within insets represent the fitted Drude profiles. Lower left panels 
show the 2D spectra of target galaxies produced by IRC data reduction pipeline. The 2D spectra have the wavelength range of 2 - 5 /xm 
which does not match the wavelength ranges of the ID spectra presented. For these 2D spectra, the shorter wavelength end is on the right 
side opposite to the ID spectra. Lower right panels show ancillary images of either the R-, or Spitzer IRAC 3.6 /im band images. Six 
galaxies with Spitzer IRAC 3.6 /im band images are FBQS J0216-0444, 2MASX J02165778-0324592, 2MASX J02191605-0557269, 2MASX 
J02193906-0511336, 2MASX J02195305-0518236, and 2MASX J16182316+5527217. Other than these six galaxies, all galaxies are with 
R-band images. These ancillary images are aligned to the 2D spectrum images to have the same spatial direction. Within the images, the 
extraction regions are drawn to the scale. 
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Fig. 3. — The stacked spectra of subsamples. (a) The spectrum stacked with the entire sample of 20 galaxies, (b) The spectrum spectrum 
stacked with the seven starburst SED galaxies. The normalized stacked spectra are overplotted for the entire sample and the starburst 
SED subsample respectively. These normalized stacked spectra are the mean of normalized individual spectra which are normalized by 
continuum flux within the wavelength range of 3.0 ~ 3.6 fim in the rest frame. We adjust the flux levels of these normalized stacked spectra 
in order to avoid overlap between the spectra. We also show the fltting results for the stacked spectra. Overplotted blue and redlines, and 
the inset are the same with Figure [21 
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Fig. 3. — The stacked spectra of the subsamples. (c) The spectrum stacked with the 12 AGN SED galaxies, (d) The spectrum stacked 
with the 15 non-detection galaxies. Overplotted blue spectra represent the normalized stacked spectra. 
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Fig. 4. — Correlation between the luminosity of the 3.3 /im PAH line and infrared luminosity. Bigger filled circles represent our sample. 
Blue circles represent galaxies with starburst SEDs while red circles represent galaxies with AGN SEDs. A green circle represent a galaxy 
with composite SE P. T he symbols without error bars, rather with leftward arrows represent upper limits. We describe how to decide 
these upper limits in l4.1l Overplotted are the data from Imanishi et al. ( 2008), Saiina et al. ( 2009), Rodriguez-Ardila &; Viegas (2003), and 
Umanishi et al. (2010). The subsamples of Imanishi et al. (2008) and Imanishi et al. (2010) are represented by small empty circles and small 
empty squares, respectively. Their subsamples with AGN SEDs are represented by red symbols and their subsamples with Hll-like SEDs 
are represented by blue symbols. Green symbols represent their subsamples with composite SEDs, or of uncertain SEDs. For the samples 
of llmanishi et a l. (2008) and Imanishi et al. (2010), data points with leftward arrows are upper limits as well. Small red crosses represent 
the sample of [Rodriguez-Ardila &; Viegas (2003), while small red empty triangles represent the high redshift sample of Saiina et al. (200^). 
The black solid line shows the linear fit for the 3.3 //m PAH emission detection from the combined sample of AMUSES and literature. The 
fit gives the the correlation between LpAH3.3iJ,m and Ljr as log (Lir) = (1.16 ±0.30) x log (LpAHS.s) - (3.11 ± 0.34). Then, the dotted 
line shows the linear fit to the data points with SB SEDs. The fit gives the correlation as log (Lir) = (1.37 ±0.17) x log (LpAHS.a) - (12.18 
± 0.75). Also presented is the fit to the non-ULIRG sample with SB SEDs. This fit gives the correlation as log (i^m) = (0.82 ±0.13) X log 
(LpAHS.s) ± (10.58 ± 0.86). It is shown as the dashed line. 
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Fig. 5. — Correlation between luminosities of the 6.2 fim PAH emission feature and infrared luminosity . The symbols are same with 
those of Figure |4] Again, the solid line is the linear fit which gives the correlation between these luminosities for the AMUSES sample. It 
gives log (Lir) = (0.95 dz 0.03) x log (Lpah6.2) + (4.21 ± 0.19), where Ljr and Lpah6.2 are in the unit of erg sec"-*^. 
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Fig. 6. — Correlation between luminosities of the 3.3 /im PAH emission feature and the 6.2 fim PAH emission feature. The symbols are 
same with those of Figure [H Again, the black line is the linear fit which gives the correlation between these two PAH luminosities for all 
the data points excluding the upper limits. It gives log (Lpah6.2) = (0.83 ± 0.06) x log (Lpahs.s) + (7.88 ± 0.41), where LpAH3.3 and 
-^PAH6.2 are in the unit of erg sec"-*^ The small inset shows 66 sources with L/i^ > 10-*^^ Lq. The fit within the inset is same with the one 
given for the entire sample. 



